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ABSTRACT 

Context. High angular and spectral resolution observations can provide us with fundamental clues to the complex circumstellar 
structure of young stellar objects (YSOs) and to the physical processes taking place close to these sources. 

Aims. We analyse the main physical parameters and the circumstellar environment of the young Herbig Be star HD 98922. 

Methods. We present AMBERA^LTI high spectral resolution (R=12000) interferometric observations across the Bry line, accom¬ 
panied by UVES high-resolution spectroscopy and SlNFONl-AO assisted near-infrared (NIR) integral field spectroscopic data. To 
interpret our observations, we develop a magneto-centrifugally driven disc-wind model. 

Results. Our analysis of the UVES spectrum shows that HD 98922 is a young (~-5xl0^ yr) Herbig Be star (SpT=B9V), located at 
a distance of 440±jg pc, with a mass accretion rate (M„cc) of ~(9±3)x 10“’ Mo yr“*. SINFONI K-band AO-assisted imaging shows a 
spatially resolved circumstellar disc-like region (~140 AU in diameter) with asymmetric brightness distribution. Our AMBERA’LTI 
UT observations indicate that the Bry emitting region (ring-fit radius ~0.31±0.04 AU) is smaller than the continuum emitting region 
(inner dust radius ~0.7±0.2AU), showing significant non-zero V-shaped differential phases (i.e. non S-shaped, as expected for a 
rotating disc). The value of the continuum-corrected pure Bry line visibility at the longest baseline (89 m) is ~0.8±0.1, i.e. the Bry 
emitting region is partially resolved. Our modelling suggests that the observed Bry line-emitting region mainly originates from a disc 
wind with a half opening angle of 30°, and with a mass-loss rate (M„.) of ~2xl0“’ Mq yr“*. The observed V-shaped differential phases 
are reliably reproduced by combining a simple asymmetric continuum disc model with our Bry disc-wind model. 

Conclusions. In conclusion, the Bry emission of HD 98922 can be modelled with a disc wind that is able to approximately reproduce 
all interferometric observations if we assume that the intensity distribution of the dust continuum disc is asymmetric. 

Key words, stars: formation - stars:circumstellar matter - stars: pre-main sequence - stars: variables: Herbig Ae/Be: individual 
objects: HD 98922 - techniques: interferometric - techniques: high angular resolution 
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* Based on observations collected at the VLT (ESO Paranal, Chile) 
with programmes 075.C-0637(A), 083.C-0236(A-D), 090.C-0192(A), 
090.C-0378(A) and 090.C-0371(A) 


1. Introduction 

Herbig Ae/Be stars ar e intermediate-mass pre-main-sequenc e 
stars (2-10 Mq; see e.g. lHerbigll960tlWaters & Waelkensll998h . 
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They are the intermediate-mass counterparts of Classical T Tauri 
stars (CTTSs) and are a fundamental link between low- and high- 
mass star formation. As is the case for their low-mass counter¬ 
parts, the circumstellar environment of Herbig Ae/Be stars is 
not fully understood. In particular, their innermost disc regions, 
within 1 AU from the central sources, harbour both accretion and 
ejection processes, whose study is of fundamental relevance to 
understanding how stars form. Herbig Ae/Be stars display sev¬ 
eral features indicative of their complex circumstellar environ¬ 
ments and of the accretion and ejection processes taking place 
close to the sources: infrared excess ind icating circumstellar 
mater ial and disc, UV excess and veiling (iDonehew & Britt^ 
I 2 OI ih indicative of accretion, as well as several permitted and 
forbidden lines which can be used to trace both accretion and 
ejection processes (e.g. Hi, Hei, or [Oi], [Feii] lines). 

However, owing to the small spatial scales involved (1 AU 
corresponds to 10 milliarseconds -mas- at 100 pc) high angular 
resolution is required to spatially resolve the inner disc region. 
Additionally, to study the gas kinematics and disentangle the dif¬ 
ferent processes, high spectral resolution is needed as well. Near- 
IR interferometry at medium and high spectral resolution has 


processes in the inner disc re 

gions (see e. 

2 . lEisner et al. 2007t 

IKraus et al.ll2008l:lEisner et al .1120091 20141 

Weigeltetal.ll201D. 


Observations have been mostly focused on the Bry line 
(2.166/im), which is bright in Herbig Ae/Be stars and can 
often be associated with accretion or ejection processes. No¬ 
tably, only five Herbig Ae/Be stars have been observed so far 
thro ugh high spec t ral resolution NIR interfe r ometry (R ~ 12000 ; 


see IWeigelt et al.l 1201 It 

Kraus el alJ 20121 

iGarcia Lonez et al.l 2015 

; lEllerbroek et al.l 


Garcia et ^ 1201 .3^ 


tions very well (see Grinin & " 

ambovtsevall201 ItjWeiffelt et al.l 

1201 UlTambovtseva et al.l2014l; 

Garcia Lonez et al.l2015). In this 


from a stellar wind, X-wind, or magnetospheric accretion, and 
they favour the last scenario. Moreover, a strong wind/outflow 
activity is also observed fro m the P Cygni profi les of the Ho- and 
Sin lines {Vr„d ~ 300kms~*: lGradv et al.ll 19961) at optical wave¬ 
length s, as well as from the Hei line profile (lOudmaiier et al.l 
I 2 OI Ih in the near-IR. 

Visual extinction values towards the object are between 
0.3 and 0.5 mag, a nd the average mag nitude in the V band is 
6.67 mag (see e.g. iManoi et al.l l2006h . The inclination of the 
HD 98922 system axis with respect to the plane of the sky is 
poorly known. A ver y rough estimate of 45° was derived in 
Blondel^Djie ( 2006h by modelling the UV spectrum, whereas 


lar, VLTI/AMBER observations of the Bry at high spectral 
resolution have been proved effective in resolving the disc- 
wind region and in inferring the flux contribution of the un¬ 
resolved magnetosphere or X-wind region to the total line 
flux. Our group has successfully developed disc-wind, X- 
wind, and magnetospheric models that match these observa- 


paper we study the circumstellar environment of HD 98922, 
through VLTI/AMBER interferometry, along with UVES (high- 
resolution UVB spectroscopy) and SINEONI-AO assisted (NIR 

integral field spectroscopy) ancillary data. _ 

HD 98922 is a Herbig Be star of spectral-type B9V (iHoukI 
[1921), although a more recent analysis su ggests it might be a 
later spectral type ('A2III: lHales et al]l2014l) . Its SED has a large 
infra red excess, indicati ng the presence of an extended dusty 
disc (iMalfait et alJll998h . Several optical and near-IR emission 
lines prove its strong circumstellar activity. In particula r , from 
the analysis of the [Oi] emission line profiles- fAcke et ari(l2005h 
conclude that this line cannot originate from the self-shadowed 
dusty disc surface, but likely from a rotating gaseous disc in¬ 
side the_du^t;sublimaton£^us. Erom the analysis of the Bry 
line. lGarcia Lopez et alJ(l2006l) derive a high mass-accretion rate 
{Mace ~ 2 X 10“^ Mq yr“\ assuming a distance of 540pc), indi¬ 
cating that the Herbig Be star is still actively accreting. AMBER - 
LR and MR interferometric studies from iKraus et ^ (l2008h 
show that the Bry line visibility increases with respect to the 
continuum visibility. In particular, the Bry line-emitting region 
is not resolved with baseline lengths up to ~60m (an upper 


Hales et alT i 2014t) report an estimate of 20° from CO modelling. 

Einally, the distance to the source is quite uncertain. 
Trig onometric distance measure d with the HIPPARCOS satel¬ 
lite (Ivan den Ancker et alJ 1199^ gives a value of d >540 pc, 
which is well beyond the HIPPARCOS observational range 
(d ~ 350^00pc). A later revisio n of these data pr ovided 
a larger value of d - 1190/®^° (Ivan LeeuwenI 1200^ . This 
value translates into an extremely large stellar luminosity for 
HD 98922 of ~10^Lo. As a consequence its position in 
the Hertzsprung-Russell diagram is well outside the locus of 
the pre-main seque nce Herbig Ae/Be stars (see e.g. Pig. 4 in 
lAlecian et al.l l20 1 31 see also a more detailed discussion in 
Sect. 13.1.1b . namely its position cannot be reproduced by us¬ 
ing the evoluti onary tracks of He rbig Ae/Be stars. Recently, 
the analysis of iHales et al.l (l2014l) has provided a closer dis¬ 
tance of 507 /;[q 4PC. It is also worth noting that this sour ce is 
a wide binary (~8", PA - 343° iDommanget & N^ll994l) and 
it is possible that a second closer companion might also be 
presen t, as hinted by spe ctro-astrometric observations (<0'.'5, 
PA~0° iBaines et al ]|200a . 

The paper is structured as follows. Spectroscopic and in¬ 
terferometric observations and data reduction are presented in 
Sect. 12] Section |3] reports results from UVES high-resolution 
spectroscopy, along with the derived stellar physical parame¬ 
ters and distance. Section|4|presents the results of our SINPONI 
observations. We present our interferometric results across the 
Bry line, along with its geometrical and physical modelling in 
Sects. |5] and |6| Section Q provides a discussion on the origin 
of the Br y emitting region and of the asymmetries observed in 
the interferometric observables. Einally, our conclusions are pre¬ 
sented in Sect.|8] 


2. Observations and data reduction 

2.1. VLT/UVES spectroscopy 

HD 9 8922 was observ ed on 21 March 2005 with 
UVES (iDekker et al.ll2000h . the echelle spectrograph mounted 
on the VLT/UT2. The observations were performed with the 
blue arm, covering a spectral range between ~3750A and 
~4990 A. A slit width of 074, which gives a spectral resolution 
of 80000, was adopted, and the total integration time was 
100s. Wavelength calibration frames were taken with a long slit 
and a ThAr arc lamp. The data were reduced using the UVES 
pipeline v3.2.0 (iBallester et al.ll2()()0h available from the ESO 
Common Pipeline Library. We converted the wavelength scale 
to the heliocentric rest frame. No standard star is available to 


limit of 1 mas, or 0.5 AU at a distance of 540 pc, is reported). 

provide an accurate flux calibration. Owing to the ver 

V small 

whereas the continuum-emitting region ring-fit diameter is es- 

photometric variabilitv of HD 98922 (Ide Winter et al. 

120011). 

timated to be 4.6+0.1 mas (or 2.48+0.05 AU at a distance of 

we adopt the photometric value reported in Manoi et al 

(2006h 

540pc). The authors conclude that the Bry line could originate 

(otb= 6.81 mag) to calibrate the spectrum. 
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2.2. SINFONI AO spectral-imaging 

HD 98922 was observed on 2 February 2013 with the 
VLT/ SINFONI integral field spectrograph (lEisenhauer et al.l 
12001 in the K-band (1.95-2.45//m) at the highest spectral 
(R ~4000) and spatial resolution (12.5x25 mas pixel scale). The 
detector integration time (DIT) and number of sub-exposures per 
frame were 2 s and five, respectively, both for the on-source (tar¬ 
get) and the off-source (sky) position. In total, eight target-sky 
pairs were obtained at position angles (PAs) 0°, 90°, 180°, and 
270° (two pairs per PA), corresponding to a total integration time 
on-source of 80s. These observations were AO-assisted (using 
HD 98922 as a natural guide star), covering a field of view (FoV) 
of0.8"x0.8". 

To correct for atmospheric transmission, observations of a 
telluric standard star of spectral type A were also performed. 
The main data reduction process was d one using the SINFON I 
data-reduction pipeline in GASGANO (iModigliani et al.ir2007h . 
i.e. dark and bad pixel masks, flat-held corrections, optical dis¬ 
tortion correction, and wavelength calibration using arc lamps. 
As a result, we obtained four calibrated datacubes, one for each 
observed PA. 

We then measured the full width at half maximum (FWHM) 
of both target and standard star on the final datacubes to esti¬ 
mate the achieved spatial resolution, which is ~60mas. From 
each datacube the stellar spectrum was then extracted, and con¬ 
tinuum and Bry line images were created by collapsing the dat¬ 
acubes along the z-axis, from 2.15 to 2.16//m and from 2.164 to 
2.168 pm, respectively. 

2.3. VLTI/AMBER/FINITO Interferometry 

The log of our interferometric observations is rep orted in Ta¬ 
ble!^ HD 98922 was observed with AMBER (iPetrov et al.l 
l2007h . the NIR beam-combiner of the Very Large Telescope In¬ 
terferometer (VLTI) operated by ESO, during two different runs 
on 26 December 2012 and 28 Eebruary 2013. On both occa¬ 
sions we employed AMBER’s high spectral resolution mode in 
the /T-band (HR mode; R = 12 000 or Av ~ 25 km s ') cover¬ 
ing the spectral range from 2.147 to 2.194//m, centred on the 
Bry line emission (at 2.166jum). Observations were conducted 
with the UT2-UT3-UT4 telescopes, using the fringe tracker 
EINITO dGai et alJl2004l) for co-phasing and a detector integra¬ 
tion time of 1.0 s per interferogram. In the first run 100 inter- 
ferograms were recorded, whereas ~1500 were recorded in our 
second run. Stars HD 60276 and HD 103125 were observed with 
the same observational settings and were used as interferometric 
calibrators to derive the transfer function for the first and sec¬ 
ond run, respectively. The calibration datasets for the first and 
second runs consist of 200 and 980 calibrator interferograms, 
respectively. 

To reduce our interferograms, we used our o wn data re¬ 
ductio n software based on the P2VM algorithm (iTatulli et alJ 
l2007h . which provides us with wavelength-dependent visibili¬ 
ties, wavelength-differential phases, closure phases, and wave¬ 
length calibrated spectra. 

The fringe tracking performance of EINITO typically differs 
during the observations of target and calibrator. Therefore, we 
used the archival low-resolution visibilities of HD 98922 (216 
different visibility uv points, mostly taken during the AMBER’s 
GTO) to calibrate the continuum visibilities at the six iiv points 
of our HR observations. This calibration method provides us 
with a visibility error of about +5% due to the large number of 
employed LR visibility measurements. By using low-resolution 



spectroscopic data taken at different epochs to calibrate our data, 
we are implicitly assuming that both size and morphology of the 
object have not changed more than the data uncertainties allow. 

The wavelength calibration of the AMBER data was refined 
using the n umerous telluric lin es present in the 2.15-2.19;um 
region (see IWeigelt et al.l 1201 iL for more details on the wave¬ 
length calibration method). We estimate an uncertainty in the 
wavelength calibration of ~0.2 A or ~3 km s *. 

3. Results from LIVES spectroscopy 

3.1. Line recognition 

The UVES flux-calibrated spectrum of HD 98922 is presented 
in Eigure [T] and shows the descending stellar continuum with 
several strong absorption lines and few faint emission features. 
The same spectrum, normalised to the continuum, is shown in 
Eigure |2] (black line). The most prominent features, both in ab¬ 
sorption and emission, are also labelled in Eig. |2] 

The most prominent absorption features are the Hi lines 
from the Balmer series (from H 11 to H/l), which partially show 
emission due to the circumstellar environment. These lines show 
P Cygni profiles (see e.g. the H/3 line vs. its theoretical profile in 
Eig. |2] the peak in emission is blue-shifted with respect to the 
zero velocity, and the blue-shifted absorption is deeper than the 
theoretical value), indicative of winds or outflows. The remain¬ 
ing absorption features are mostly photospheric from atomic 
lines like Hei, Can, Eei, Eeii, Mgi, Mgii, Sii, and Sin. Other 
emission features, mostly Ee ii and Ti ii lines, located between 
~4200 and 4600 A, arise from circumstellar activity. These lines 
are characterised by a double-peaked profile, are likely pro¬ 
duced by self-absorption of the photospheric feature, and have 
a broad pedestal (full width zero intensity, FWZI ~4^.4 A, i.e. 
Av ~270-290kms“'), which originates in a region with a large 
radial velocity structure, likely located in the magnetospheric ac- 
creti on funnels or in the hot inner region o f the circumstellar 
disc (iBeristain et al.lll99^lAspin et al.l 20101) . The strongest Ee ii 
emission line (at 4925 A) shows a P Cygni profile. 

Prom the analysis of the H 11 line, we derive a value of — 
8kms“' for the stellar radial velocity (Vrad) with respect to the 
local standard of rest. Notably, the observed spectral lines do not 
show any indication of binarity: no splitting and no additional 
absorption lines are detected in the spectrum. 

3.1.1. Stellar physical parameters and distance 

To allow both identification and analysis of the observed 
line profiles we used synthetic spectra from the PHOENIX 
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Table 1. Log of the VLTI/AMBER/FINITO observations of HD 98922 and calibrators. 


HD 98922 

Time [UT] 

Unit Telescope 

Spectral 

Wavelength 

DIT" 

N" 

Seeing 

Baselines 

PAs 

Obseiwation 

Date 

Start 

End 

array 

mode‘s 

range 

(^m) 

(s) 


(") 

(m) 

(°) 

2012 Dec 26 

09:02 

09:20 

UT2-UT3-UT4 

HR-K-F 

2.147-2.194 

1 

100 

0.7-0.S 

44.1/60.6/89.4 

219.3/103/76.7 

2013 Feb 28 

04:50 

05:32 

UT2-UT3-UT4 

HR-K-F 

2.147-2.194 

1 

1540 

0.9-1.2 

43.9/61/89.4 

220.9/104.9/78.8 


Calibrator 

Name 

Date 

Time [UT] 

Start End 

UT array 

Spectral 

mode‘ 

Wavelength 

range 

Oim) 

DIT^ 

(s) 

N" 

Seeing 

(") 

Uniform-disc 

diameter^ 

(mas) 

HD 60276 
HD 103125 

2012 Dec 26 

2013 Feb 28 

07:50 

04:08 

08:03 

04:42 

UT2-UT3-UT4 

UT2-UT3-UT4 

HR-K-F 

HR-K-F 

2.147-2.194 

2.147-2.194 

1 

1 

200 

70 

0.9-1.0 

0.9-1.5 

0.92±0.06 

0.88±0.06 


Notes. Detector integration time per interferogram. Number of interferograms. High spectral resolution m ode in the Ai-band using the 
fringe tracker FINITO. UD diameter taken from JMMC Stellar Diameters Catalogue - JSDC dLafrasse et al.ll20l5) . 


library (iHusser et al.l 1201 3h . which provides high-resolution 
(R=500000) synthetic spectra based on the stellar atmosphere 
code PHOENIX. The spectral resolution was then reduced to 
that of the UVES observations by convolving the synthetic spec¬ 
tra with a Gaussian profile (IRAE task gauss). We constructed 
a grid of spectra and varied different stellar parameters, namely 
Teff (from 9000 to 11 000K, AT - 200K), logg (from 3.0 to 
4.5, Alogg = 0.5), and [Ee/H] (from -1 to 1, A[Ee/H] = 0.5). Both 
observed and synthetic spectra were then normalised to the con¬ 
tinuum. Our own IDE programme was used to match the UVES 
spectrum with the best-fitting model in three steps. Eirst, by com¬ 
paring the wings of the Balmer lines the stellar gravity was in¬ 
ferred. In this case, Balmer lines from HI 1 to Hy were used (the 
HjS line was excluded because it is too affected by circumstel- 
lar emission). Second, the effective temperature of the star 
was obtained by comparing observed and theoretical equivalent 
widths (EWs) of lines sensitive to T^s, namely the Hei lines (at 
4026, 4120, 4145, 4470 A), and the Tin line at 4176 A. Einally, 
the metal abundance was inferred by comparing the iron abun¬ 
dances for a well-dehned set of Ee i and Ee ii lines. As a result, 
we obtained the stellar effective temperature, gravity, and metal- 
licity, namely Egfr = 10400+200K (i.e. spectral type B9V), logg 
- 3.5+0.2, and [Fe/H] - -0.5+0.2. The best-fitting synthetic 
spectrum, normalised to the stellar continuum, is shown in Eig- 
ure |2] (red solid line) superimposed on the observed spectrum 
(in black). Eor comparison, Eig.|2]also shows a normalised syn¬ 
thetic spec trum (green dashed line) with the stellar parameters 
derived bv iHales et al.l (l2014l) . namely Egfr = 9000 K (SpT=A2 
III), logg^3.0, and [Fe/H]^-0.5. 

To derive th e rema ining stellar parameters, we follow 
iMontesinos et al.l (l2009h . We place the estimated values of 
the stellar parameters on a log Teff, log g HR diagram, us- 
ing evolutionary tracks and isochrones from ISiess et ^ (l2000h 
with [Fe/H]--Q.5. This provides us with an estimate of 
the stellar mass (M.=5.2+0.2Mo) and age (~5xl0^yr) for 
HD 98922, and, in turn, stellar luminosity (L, = 640+[ggLo), 
radius (R, = 7.6+° ® R©), and visual absolute magnitude (Mv=- 
1.75+0.3 mag) can be inferred from a (logTeff ''■s- logL,) HR 
diagram. 

Eollowing lMontesinos et all (l2009h . we also aim at deriving 
a better estimate of the distance (d) to HD 98922, from the dis¬ 
tance module equation: 


d — IQ(>’'^v-Av-Mv+5)/5 
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Table 2. HD 98922 stellar and disc parameters 


Stellar Parameter 

Value 

Reference 

Distance 

440+“ pc 

1 

Teff 

10400+200K 

1,2 

S pType 

B9V 

1,2 

logg 

3.5+0.2 

1 

[Fe/H] 

-0.5+0.2 

1 

M. 

5.2+0.2M© 

1 

R , 

7.6+°° R© 

1 


6:40+130 T 

1 

age 

5x10° yr 

1 

My 

-1.75+0.3 mag 

1 

my 

6.77 mag 

3 

Av 

0.3 

3 

niK 

4.28 

4 

Fobs(K)/F,(K) 

4.3 

1 

i 

20° 

5 

^rot sin i 

50+3kms-i 

6 

^acc 

(9+3)xl0-’M©yr-* 

4+ 1 


Refer ences. (1) This work; (2) iThe^alJ (Il994li: ('31 iMmoi et all 
l l2Ci03): (41 iGarcia Lopez et al.l ( l2006h : (51 IHales et al.l ( 120141) : 
(bllAlecian et^ jMir 

We a ssume mv=6.77mag and Ay=0.3 (iGarcia Lonez et all 
l2006h . and we obtain ^f=440+gQpc. By adopting this distance 
value, an estimate of the mass accretion rate {Macc) can be also 
inferred from the EW of the Bry line after correcting for the 
vis ual extinction and self-abs orption from the photosphere (see 
e.g. IGarcia Lonez et al.ll200^ . By using their published data we 
obtain Mace~(9+3)x 10“’ M© yr“*. 

All stellar and disc parameters used in this paper, derived in 
this work or adopted from the literature, are listed in Table |2 


4. Results from SINFONI-AO spectral imaging 

Our SINEONI K-band spectral images resolve a bright asym¬ 
metric scattered emission (diameter ~320 mas at 3cr, or ~ 140 AU 
at a distance of 440 pc) that extends beyond the symmetric point 
spread function (PSE; EWHM=60 mas) of HD 98922 (see the 
upper panel of Eigure|3]and its contour levels). This structure is 
not an artefact and because it is visible in all the spectral images 
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47(X) 4750 4800 4850 49(X) 4950 5000 
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Fig. 2. UVES/VLT high-resolution spectrum of HD 98922 normalised to the continuum (black solid line). The red solid line shows the best¬ 
fitting synthetic spectrum, SpT=B9\, logg=3.5, [Fe/H]=-0.5. Eor comparison, a synthetic spectrum with parameters S p7’=A2III, logg=3.0, 
[Fe/H]=-0.5 (green dashed line) is also shown. The most prominent spectral features are labelled. 


and is co-rotating in each image FoV according to the different 
observed PAs. 

The bright circumstellar structure most likely originates from 
the scattered light of the disc. Of particular interest is its asym¬ 
metric morphology. The brightness distribution is roughly repre¬ 
sented by a symmetric PSF plus a bright elliptical arc towards the 
north-north-east (NNE) (i.e. a horseshoe-like region), extending 
from PA~290° to ~110° (see upper panel of Fig. [3]l and located 
at ~90 mas (or ~40 AU) from the source. This region is about 2.5 
times brighter than its symmetric counterpart located towards 
south-south-west (SSW). The value of this ratio was obtained 
after integrating the flux in two regions (positioned outside the 


PSF, i.e. 60 mas), which have the same area and are symmetri¬ 
cal with respect to the central PSF. It is also worth noting that 
the brightness distribution of the NNE arc is not homogeneous; 
it shows an elongated peak between PA~23° and PA~65°, and 
it is up to three times brighter than the whole arc-shaped fea¬ 
ture. The total flux of this area is about one tenth of the flux 
of the central PSF Because of its elongated arc-shape, which 
cannot be fitted with the instrumental PSF, we exclude that this 
structure can be the unresolved PSF of a possible companion. 
Remarkably, our SINFONI-AO continuum images do not show 
any companion between ~800 mas (or ~350 AU) and ~60 mas 
(or 26.4 AU) from the central source. The lower limit of the flux 
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Fig. 3. Upper Panel: SINFONI K-band continuum image (between 2.15 
and 2.16 pm). Labelled fluxes are in W m“^ pm arcsec^^. Contour levels 
are 4, 5, 50, 150, 250, 400, 500, and 1000 cr. Lower Panel: Intensity 
ratio image of two SINFONI K-band continuum images (between 2.15 
and 2.16^m) taken at PA of 0° and 180°, showing the asymmetric disc 
emission. Contours show the K-hand continuum of the PA=0° image 
(see upper panel), and give the position of both symmetric PSF in the 
centre, and asymmetric disc emission. Contour levels are as in the con¬ 
tinuum image. 


density in our data is ~5xl0“*' ergs^' cmT ^which cor¬ 
responds to K ~ 9.8 mag or Mk ~ 1.6 mag at the estimated 
distance of 440pc. Adopting the lSiess et al.l (l2000l) evolutionary 
tracks and assuming an age of ~5xl0^yr for any possible un¬ 
detected companion, the Mk estimate provides us with an upper 


limit constraint on the spectral type and mass of SpT~M0 and 
M.~0.5Mo. 

To better display the observed brightness asymmetry in the 
disc, the lower panel of Figure [3] shows an intensity ratio image 
that is computed from two SINFONI K-band continuum images 
(between 2.15 and 2.16/rm) with PAs of 0° and 180°, respec¬ 
tively. The central unresolved PSF disappears (ratio value ~1), 
whereas the bright NNE arc becomes clearly visible (ratio values 
between 1.5 and 5) in the image. Apart from the aforementioned 
elongated peak, there is a dip at PA~0° and a second arc-shaped 
peak towards the north-west that is less bright. 

It is worth mentioning that the PIONIER/VLTI reconstructed 
image s of the HD 98922 inner disc (see Eigure 2 in iKluska et al.l 
l2014l) . show a similar brightness disc asymmetry, namely the 
north-eastern side of the disc is brighter than the south-western 
side. These interferometric data were acquired between 20 De¬ 
cember 2012 and 20 Eebruary 2013 (Kluska priv. comm.), si¬ 
multaneous with our AMBER and SINEONI observations. The 
observed PIONIER disc has, however, a size of ~3-4mas (~1- 
1.5 AU), which is well within the (spatially unresolved) PSE. 
This suggests that the disc might have a similar morphology 
from ~1 to 60-70 AU from the source and that the origin of this 
brightness asymmetry might be the same. 

The Bry line is the only emission line detected over the K- 
band continuum of the analysed SINEONI spectral-images. This 
emission is not spatially resolved, being enclosed within 60 mas 
(or ~26 AU) from the central source. 


5. Results from VLTI/AMBER interferometry and 
geometric modelling 

5.1. Interferometric observables: Visibilities, differential 
phases, and closure phase 

Our VLTI/AMBER observations provide us with four direct 
observables, namely the Bry line profile, visibilities, differen¬ 
tial phases, and closure phases. These observables allow us 
to retrieve information about the inner-disc emitting region of 
HD 98922, in particular the size and kinematics of the Bry emit¬ 
ting region, as well as its displacement with respect to the con¬ 
tinuum emission. 

Eirst of all, we stress that observations in both VLTI runs 
were taken with the same AMBER settings and UT configura¬ 
tions, which result in similar projected baselines and position 
angles (see Cols. 10 and 11 of Table[T]i. The results from both 
runs are identical inside the error bars. However, those from the 
second run (28 Eebruary 2013) are clearly much less noisy than 
those from the first run (26 December 2012) owing to the longer 
total integration time. Eor this reason, in this section we only 
present the results from the second run (EigurelUi, whereas those 
from the first run are shown in Appendix lAl (Eigure lA. lb . The 
different panels in Eiguresl4l and lA.il show (from top to bottom) 
the Bry line profile (upper panels), visibilities (middle upper 
panels), differential phases (middle lower panels), and closure 
phases (lower panels). 

The wavelength-dependent visibilities (middle upper panels 
of Eig.l4land lA.lb clearly increase across the Bry line-profile at 
all baselines. This indicates that the Bry emitting region is more 
compact than the continuum emitting region. Previous spectro- 
interferometric studies of HD 98922 at medium resolution (R = 
1500) ha ve also detected an increase in the visibility within the 
Bry line (iKraus et al.l l2008l) with baseline lengths up to 60 m. 
Our HR mode AMBER observations allow us, for the first time, 
to measure the visibilities and phases in ~30 different spectral 
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velocity (LSR) [km/s] 
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X [|rm] 

Fig. 4. AMBER observation of HD 98922 at spectral resolution of 
R=12 000 (28 Eeb 2013). Erom top to bottom: wavelength dependence 
of flux, visibilities, wavelength-differential phases, and closure phase 
observed at three different projected baselines (see labels in figure). Eor 
clarity, the differential phases of the first and last baselines are shifted 
by -1-30° and -30°, respectively. 


channels across the Bry line with baseline lengths up to 90 m, 
providing a higher spatial and spectral resolution. Furthermore, 
our data indicate that the visibilities from the three baselines are 
slightly asymmetric with respect to the Bry peak, namely the 


visibilities of the blue-shifted wings are systematically smaller 
than those of the red-shifted wings. This might indicate that, at 
high velocities, the red-shifted Bry emitting region is slighlty 
more compact than the blue-shifted one, or that there is a highly 
asymmetric continuum. 

Differential phases (d); lower middle panels of Fig. |4] and 
lA.ll) measure the displacement of the line-emitting region with 
respect to the continuum-emitting region. For example, the ac¬ 
curacy reachable with a 90 m baseline and with O ~10° is 
about 0.14 mas, or 0.06 AU at a distance of 440 pc. Therefore 
the observed displacements provide us with a very sensitive 
measurement of the gas kinematics on scales of a few stellar 
radii. For the first time in this object, we detect a clear pho¬ 
tocentre shift of the Bry line with respect to its continuum at 
the 61m and 89 m baselines. The line displacement is observed 
in different velocity channels (both blue- and red-shifted). No¬ 
tably, differential phases in both blue- and red-shifted channels 
have the same sign, and therefore they do not exhibit the typ¬ 
ical “S” shape usually observed in rotating discs. Additionally, 
the red-shifted channels in both baselines show larger differen¬ 
tial phases with respect to the blue-shifted ones. Because these 
are line-to-continuum phases, the observed shifts in the differ¬ 
ential phases could be caused by continuum shifts, line shifts, 
or both. The maximum value of the line/continuum displace¬ 
ment is 0.02+0.04 mas, 0.14+0.04 mas, and 0.16+0.03 mas for 
the 43.9 m, 61 m, and 89.4m baselines, respectively. 

Finally, it is also worth noting that the closure phases do not 
show any significant deviation from zero within the error bars 
(see lower panels of Figs. l^land lATI) : we measure an average 
value of -0.44°+2.55°. 

5.2. Geometric modelling: size of Bry and continuum 
emitting region 

Our AMBER high spectral resolution data also provide us with 
a measurement of the visibilities for the pure Bry line-emitting 
region in the different spectral channels across the Bry emission 
line (see Figure|5]l. The continuum-corrected line visibilities are 
fundamental for measuring the size of the Bry line-emitting re¬ 
gion. Within the wavelength region of the Bry line emission, the 
measured visibility has two components: the pure line-emitting 
part and the continuum-emitting part, which includes continuum 
emission from both the circumstellar env ironment and the un - 
resolved central star. Therefore, following lWeigelt et al.l (l2007h . 
the emission line visibility Vety in each spectral channel can be 
written as 

FBryVBry = V^tolVtotP + \F,V,\^ - 2 F,V, ■ COiO, (2) 

where Fsry is the wavelength-dependent line flux, Vtot (Etot) is 
the measured total visibility (flux) in the Bry line, Vc (Fc) is 
the visibility (flux) in the continuum, and O is the measured 
wavelength-differential phase within the Bry line. 

In our analysis, we also included the intrinsic Bry photo- 
spheric absorption feature of HD 98922, using a synthetic spec¬ 
trum with the spectral type and surface gravity values as in Ta¬ 
ble E] 

The size of the continuum was obtained using the AMBER 
GTO archival data at low resolution. By fitting a circular sym¬ 
metric Gaussian model, we derived a Gaussian half width at half 
maximum (HWHM) radius of 1.8+0.3mas (or ~0.79+0.13AU 
at a distance of 440pc). A similar value of 1.6+0.3mas (or 
~0.70+0.13 AU) was inferred by fitting a ring model with a ring 
width of 20% of the inner ring radius. 
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velocity (LSR) [km/s] 



2.160 2.162 2.164 2.166 2.168 2.170 2.172 


X [um] 


Fig. 5. Comparison of the observed (pink solid lines) and modelled 
(black solid lines) continuum-corrected (pure) Bry line visibilities of 
our AMBER observation of HD 98922. From top to bottom: wavelength 
dependence of flux and visibilities of the first, second, and third base¬ 
line. In each visibility panel the observed total visibilities (red, green, 
blue, as in Fig.|4), the observed pure Bry line visibilities (pink), and the 
modelled pure Bry line visibilities (black) are shown. 


To get reliable line visibility values, we derived the 
continuum-corrected (hereafter pure) line visibility only in the 
spectral regions where the line flux is higher than ~10% of the 


continuum flux. The results are shown in Figure |5] The aver¬ 
age line visibility is ~ 1+0.2 for the shortest baseline (~44m), 
~0.9±0.15 for the medium baseline (~61 m), and ~0.8+0.1 for 
the longest baseline (~89m). For this baseline, the line is spa¬ 
tially resolved even when the error bars are taken into account. 

The size of the line-emitting region was obtained by fit¬ 
ting a circular symmetric Gaussian model to the line visibili¬ 
ties presented in Fig. |5] We obtain a Gaussian HWHM radius 
of 0.65+0.09mas or ~0.29+0.04AU at a distance of 440pc. 
Similarly, by fitting a ring model with a ring width of 20% 
of the inner ring radius, an inner radius of 0.70+0.09 mas, or 
-0.31+0.04 AU, is inferred. 

6. Disc and disc-wind modeis 

To constrain the physical processes taking place in the inner¬ 
most disc regions, a proper physical mode lling is required. As 
in our previous interfero metric works (see IWeigelt et al.ll201 It 
iGarcia Lonez et al'1l2015h . to analyse the observed Bry line pro¬ 
file, line visibilities, differential phases, and closure phase we 
employed a previously developed disc-wind (DW) model (see 
Sect. 16.1b and developed a continuum disc model; both models 
are adapted to the HD 98922 stellar parameters and observations. 

6.1. Disc-wind modei 


A detailed description of our DW model, its parameters, and the 
algorithms used for the model computation can be found in our 


previous papers (seelGrinin & Tambovtsevall201 ll 

Weiffelt et al.l 

1201 ItlTambovtseva et al.ll2014 ; Garcia Lonez et al. 

201.511. In the 


following, we briefly outline the model’s main characteristics 
and define its free parameters, which are listed in Table |3] 


The developed model is a warm disc-wind model (see e.g. 
ISafierlFig^ iGarcia et al.ll200ll) . including only hydrogen atoms 
with constant temperature (-10000K). The wind is rapidly 
heated by ambipolar diffusion to a temperature of -10 000 K, 
and the wind electron-temperature (Te) in the acceleration zone 
near the disc surface is not high enough to excite the Bry line 
emission, which is emitted further out along the wind stream 
lines. 

Briefly, the disc wind is launched from a disc inner radius lo\ 
to an outer radius un (called wind footpoints), where N is the 
number of streamlines. Its half opening angle (0) is defined as 
the angle between the innermost wind streamline and the system 
axis. The local mass-loss rate per unit area on the disc surface is 
defined as m„{oS) - where y is the mass-loading parameter 
that controls the ejection efficiency. Tfierefore, the total mass- 
loss rate {M„) is 

M„ ~ ^ J" 2 TT w da>. (3) 

Finally, the yS parameter in Table 0 reg ulates the accelera¬ 
tion o f the wind along the streamlines (see iTambovtseva et all 
l2014h . Some model parameters such as the stellar parame¬ 
ters, and the disc inclination angle were derived in this work or 
were taken from the literature (see Tab. la. The explored range 
of values of M„ varies from -0.05 to 0.3 times the inferred 
mass accretion rate {Macc, see Table|2]l. The value for the coro¬ 
tation radius (/?co/El) obtained from the v™, sin i value (Table|2]l is 

* with Rcorr=(GM,/vl^y^^, where G is the gravitational constant, M, 
the stellar mass, and the stellar angular velocity. 
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Rcor-i-55Rf The estimated value of R^or is then set as a lower 
limit to the inner disc-wind launching radius value wi to allow 
for the presence of a hypothetical magnetosphere within this re¬ 
gion. 

To find the best model that reproduces our ob servations, we 
adopt the following steps (for more d etails, see also lWeigelt et alJ 
1201 It iGrinin & Tambovtsevall201ll) . First, we fit the observed 
Bry line profile and intensity. To this end, we solve the equations 
of statistical equilibrium and compute the population of the H i 
atomic levels and the ionisation degree a long each streamline 
using the n umerical codes developed by iGrinin & Mitskevk^ 
d 19901) and iTambovtseva et all (1200 ih for moving media. Sec¬ 
ond, we calculate the intensity of the line radiation in the entire 
emitting volume filled in by the disc wind and compare both 
modelled and observed line profiles. When a good agreement is 
found, we compute the brightness distribution of the disc wind 
and the disc continuum at different radial velocities, i.e. a two- 
dimensional intensity distribution map, which provides us with 
the interferometric observables. The modelling results are, in 
this way, directly compared with the observations until the best- 
fit model is found. The best model is that which best matches, 
within the error bars, all eleven interferometric observables (line 
profile shape and intensity, visibilities, pure line visibilities, dif¬ 
ferential phases, and closure phase) of our observation within the 
selected range of parameters listed in Tabled 

This latter reports the free model parameters, the explored 
range of parameters (Column 2), and the values of our best-fit 
model, called P5 (Column 3 of Tab. O. P5 model was selected 
after checking hundreds of different models configurations (see 
Column 2, Tab. [3]). 

6.2. Continuum disc modei 

To reproduce our AMBER HR interferometric observations, we 
have assumed that the HD 98922 continuum emission con¬ 
sists of three different components : the star, an inner disc , 
and an outer dusty disc (s ee e.g. iTannirkulam et al.l l2008bl: 
iDullemond & MonnieJ20f^ . 

The main constraints for our developed temperature-gradient 
model of the continuum emission are provided by: i) the 
^obs{K)l¥t{K) ratio, namely the ratio between the HD 98922 ob¬ 
served flux continuum in the /T-band (star plus accretion disc) 
and the theoretical stellar flux, and ii) the absolute visibilities 
of the continuum. The stellar continuum at the Bry rest wave¬ 
length {F,(K)) was computed from a Kurucz synthetic spectrum 
with 7’eff=10400K and logg-3.5 (see Sect. 13.1.H and Tab. |2]i, 
whereas the F„tisiK) is computed from the YSO /T-band mag¬ 
nitude (see Tab. O. The adopted disc inclination angle is the 
same as for the disc-wind model (i=20°, see Tab. H. To fit 
the disc continuum level of our observation s, we use an inner 
gaseous disc plus an outer dusty disc, as in ITannirkulam et ^ 
d2008ah and iBenistv et al.l (l2010h . The sum of fluxes from the 
gaseous and dusty disc matches the observed flux close to the 
Bry wavelength. The inner disc has constant brightness distri¬ 
bution and ranges from R,„=4R, (or 0.14AU) to Rs-2QRt (or 
0.7+0.2AU). The dust temperature (Tjo) at this radius (Rg) is 
150 0 K, that is, the silicate dust sublimation temperature (see 
e.g. IDullemond & Monnieill2010h . The dust temperature expo¬ 
nentially decreases with the distance r from the star accordingly 
with the following power law: T^ir) - TdoirjRgY, with a=-0.5. 
The intensity of the radiation in the continuum was calculated in 
the blackbody approximation. The disc brightness is distributed 
symmetrically along azimuthal direction. Fig. |6] (upper panel) 


shows the disc continuum model plus the Bry intensity distribu¬ 
tion map (for Vr=0kms ') of the best-fit disc-wind model (P5). 



X/R. 

Fig. 6. Upper panel. Bry intensity distribution map (for Vr=0kms“') 
of the best-fit disc-wind model (P5) overlapped on the map of the disc 
continuum with symmetric brightness distribution. X and Y coordinates 
are in R,. Lower panel. Same as in the upper panel, but for the disc 
continuum with asymmetric brightness distribution. 


6.3. Modeiiing resuits 

As seen in Sect. |5] the Bry emission line shows V-shaped differ¬ 
ential phases and slightly asymmetric visibilities (see Fig.|4]i. 

Because of the aforementioned shape of the differential 
phases, our first attempts at modelling these interferometric ob¬ 
servables with our disc-wind model plus the continuum disc 
model described in Sect. I6.2l were unsuccessful (see Fig. IBHI in 
Appendix |B]i. Indeed, the combination of these disc and disc- 
wind models produces, as expected, symmetric S-shaped differ¬ 
ential phases (originating from disc rotation, see Figure |7j upper 
panel) and symmetric visibilities in the Bry line. 

In principle, the unusual shapes of the observed visibilities 
and differential phases, which are a measurement of the photo¬ 
centre shifts, can be caused by an asymmetric Bry line-emitting 
region, an asymmetric continuum-emitting region, or by a com¬ 
bination of both. As a consequence, the differential phase in the 
blue or red wings of the Bry can be displaced, and the visibil¬ 
ities become asymmetric. An asymmetric Bry line-emitting re¬ 
gion could be a plausible explanation for the observed differ¬ 
ential phase. In theory, this might originate from a collimated 
wind/iet, whose red -shifted lobe is hidden by the disc (see e.g. 
IBenistv et al.ll201^ . However, in this case, we would observe a 
blue-shifted photocentre shift, whereas our observations clearly 
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Table 3. Disc-wind accretion model parameters: range of explored values and values for our best-fit model P5 


Parameter" 

Range of Values 

Model P5 

Half opening angle (0i ) 

20°-60° 

30° 

Inner radius - u)\ (R») 

2-3 (0.07-0.1 AU) 

3 (0.1 AU) 

Outer radius - w^(R*) 

4-30 (0.13-1 AU) 

30(1 AU) 

Acceleration parameter (f) 

3-7 

5 

Mass load parameter (y) 

2-5 

3 

Mass-loss rate - M„,(Mo yr^') 

10-*-3xl0-^ 

2x10-’ 


Notes. See lWeigelt et alJ (1201 ih for definition and detailed description of the listed disc-wind model parameters 


show a red-shifted peak in the differential phases. On the other 
hand, it is worth noting that the HD 98922 NIR continuum emis¬ 
sion shows an asymmetric brightness distribution on small and 
large spatial scales (from ~1 AU to tens of AUs), as observed 
in both VLTI/PIO NIER and SINFONI data (see Sect, g] and 
iKluska et al.ll20l4) . Therefore, an asymmetric continuum might 
be responsible for the observed shape in the differential phases 
and visibilities. 

To verify this hypothesis, we modified the disc model of 
Sect. 16.21 developing a very simple asymmetric continuum disc 
model, which mimics the difference in brightness between the 
two halves of the disc. However, it is beyond the scope of this 
paper to model the observed disc asymmetries in details and at 
small scales. Therefore, for the sake of simplicity, the intensity 
of the radiation was increased by a factor fc in the northern half 
of the disc and decreased by the same factor in the southern half 
of the disc of both gaseous and dusty disc regions. As a conse¬ 
quence, the disc brightness from PA 270° to 90° (i.e. the northern 
half of the disc) is brighter than the southern half of the disc (ex¬ 
tending from PA 90° to 270°). Lower panel of Figure. |6] shows 
the asymmetric disc continuum model plus the Bry intensity dis¬ 
tribution map (for v^^Okms *) of the best-fit disc-wind model 
(P5). The brightness contrast (fc) between the two halves of the 
disc is a free parameter of the model. The asymmetric disc was 
then rotated by 10° steps until we obtained a solution that best 
fits both visibilities and differential phases at the three observed 
baselines. In order to fit the continuum level of the three base¬ 
lines, we changed the fc value, which was varied from 1 to 4. 
The best continuum model has an fc value of 2.5 and a disc ro¬ 
tation of 70°. Although very simple, our asymmetric disc model, 
combined with the disc-wind model described in Sect. 16.11 is 
able to reproduce the V-shaped differential phases (see Fig. |7] 
middle panel). Middle and lower panels of Fig. Q present differ¬ 
ential phases for different rotation of the disc and the disc wind: 
counter clockwise (b; middle panel) and clockwise rotation (c; 
lower panel). It should be noted that the shapes of the differen¬ 
tial phases in the two cases are slightly different. The reason is 
the disc wind. By modifying the direction of the disc-wind rota¬ 
tion, the distribution of the disc-wind brightness changes but the 
distribution of the disc brightness remains unchanged. Thus, the 
differential phases are sensitive to the direction of the disc and 
disc-wind rotation. 

Figure [8] shows the Bry intensity distribution maps of our 
best disc-wind model (P5) in logarithmic scale, as a function 
of different Bry radial velocities. The adopted system axis in- 
clinat ion angle with res pect to the line of sight is 20° (see Ta¬ 
ble |2l|Haies^et^|20l3)- Columns of Tabled reports the values 
of the model free parameters. Figure 0 shows a comparison of 
our best model with the interferometric observables. Despite the 
simplicity of the proposed disc model, the interferometric ob¬ 


servables are approximatively well reproduced, demonstrating 
that V-shaped differential phases and slightly asymmetric line 
visibilities may originate from the uneven brightness distribu¬ 
tion of the disc. It is also worth to note that the value derived 
for fc (i.e. 2.5) is, within the error bars, identical to the aver¬ 
age brightness ratio between the northeastern and southwestern 
disc regions of the SINFONI spectral images (see Sect. HI, and 
that the PA of our asymmetric disc model (i.e. 70°) is compat¬ 
ible with the observed asymmetric disc orientation in Fig. [3] In 
Fig. |5] we compare observational and model results of the Bry 
pure line visibilities (pink and black solid lines, respectively). 
Our model approximatively reproduces the observed values at 
the shortest and medium baselines (upper and central panels), 
but it slightly overestimates the pure line visibility at the longest 
baseline (lower panel). This result represents a good compro¬ 
mise. To obtain a lower value of the modelled pure line visibility 
at the longest baseline, it would be necessary to increase the size 
of the disc-wind region (i.e., the value un)- On the other hand 
this would increase the values of the differential phase, and the 
resulting model would not match our observations. 

Moreover, we also note that the Bry emitting region is spa¬ 
tially extended. As a consequence, it cannot be modelled with 
a compact, spatially unresolved com ponent alone, such as a 
magnetosphere (size <2-3 R,,; see e.g. lTambovtseva et al.ll201^ 
iGarcia Lopez et aHl2015h . that cannot be spatially resolved at 
these baselines. Therefore we conclude that the Bry emission 
must be at least a combination of a compact (i.e. spatially un¬ 
resolved) emission (magnetosphere) plus an extended spatially 
resolved emission (i.e. the disc wind). The addition of a com¬ 
pact region (i.e. visibility value = 1 at all baselines) implies that 
the size of the spatially resolved component must be enlarged to 
match the observed visibility values. This can be achieved with a 
larger disc-wind footpoint radius which, in turn, gives larger 
values of the differential phases. In addition, changing the foot- 
point radius modifies the line profile as well. Thus our model 
has to take into account these three observables (visibility, differ¬ 
ential phase, and line profile). In conclusion, by comparing our 
model with the observations, we do not rule out that some of the 
Bry emission may originate from the magnetosphere (i.e. a spa¬ 
tially unresolved region), however we can exclude that this latter 
is t he main mechanism, whic h produces the Bry emission (see 
also lGarcia Lopez et al.ll20T5l) . 

Finally, it is worth comparing the mass-loss rate derived from 
our model with the mass accretion rate inferred from the ob¬ 
servations, as their ratio (M„IMacc) usually provides constraints 
on the efficiency of the accretion/ejection mec hanism in young 
stars, typically rang ing from 0.1 to 0.3 (see e.g. lShu et al.lll988[ 
iFerreira et al.ll200^ . In our case, the accuracy of the Macc value 
is ~30%. This uncertainty mainl y originates from the employe d 
empirical relationship (see e.g. iMendigutfa et al.ll201 iL l2013h . 
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Fig. 7. Differential phases for the three observed baselines (as in Fig. 0, 
calculated for the best disc-wind model P5 plus: symmetric disc model 
(a, upper panel), asymmetric disc model with a counter clockwise (b, 
middle panel), and clockwise (c, lower panel) disc and disc-wind rota¬ 
tion. For clarity, the differential phases of the first and last baselines are 
shifted by -1-20° and -20°, respectively. 


On the other hand, it is more difficult to estimate the accuracy 
of the M„ value as inferred from our best-fit model, but it is 
probably of the same order as the Mace value, owing to the use 
of our simplified model of the disc wind and to the uncertainties 
in the fit of the Bry line profile and intensity. As a result, their ra¬ 
tio is well within the expected range of values. It should also be 
noted that in our best-fit model of HD 98922 the disc wind starts 
in the close vicinity of the star. This can mean that the wind at 
least partially originates as a result of interaction of the accretion 
disc with the magnetosphere. In this case the MwjMacc ratio can 
be higher than that used in the theory of the magneto-centrifugal 
disc wind (0.1), for example due t o the effect of the mechanis m 
of the magnetic propeller (see e.g. llllarionov & SunvaevlfT975h . 


7. Discussion 


7.1. An asymmetric disc emission 

Asymmetries in the disc brightness distribution are not uncom¬ 
mon as they have been observed at NIR wavelengths in several 
Herbig Ae/Be stars and CTTSs with the Hubble Space Telescope 
(HST) or with ground-based telescopes thr ough AO assisted 
imaging or interferometric im aging (see e.g. iKrist et aTI I2005I: 
iDullemond & Monnied l2010h . These images of spatially re¬ 
solved discs displ ay a wide range of spectacular asymmetr ies, in¬ 
cluding arcs ( e.g. iFuka gawa et ^l2006l: iKrist et a Hl2012|). gaps 
and dips (e.g. IKrist et al.l 2005 ). warps (e.g. Golimowski et al.l 
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Fig. 8. Bry intensity distribution maps of our best disc-wind model P5 
(see values of parameters in Table[3 Column 3) in logarithmic scale (ar¬ 
bitrary units). Each panel shows the intensity map for a different radial 
velocity, indicated by a white label in units of kms“*. X and Y coor¬ 
dinates on each map are in R,. The system axis inclination angle with 
respect to the line of sight is assumed to be 20° (see Table0. Contin¬ 
uum emission from the disc and central star, located at coordinates 0,0, 
are not shown. 


200^ M^no^t aH 2()jj|k or more complex morphologi es (e.g. 


The origin of such asymmetries in discs varies from case 
to case, but mostly involves disc geometry, illumination (e.g. 
only one side of th e illuminated inner rim is observed; see 
iMonnier et al.l l2()(^ . dust grain composition and distribution, 
or tidal effects: a) brightness asymmetries between the two 
halves of the disc may originate from different emission in the 
forwards/backwards scattered light in incline d discs (see e.g. 
IKrist et ^l2005l:rDullemond & Monnied2010l) : b) gaps and dips 
in discs can be produced by shadowing ef fects caused by o pti- 
cally thick clouds of dust and gas (see e.g. IKrist et al.ll2()()5h : c) 
arcs of different brightness and size may originate from shadow¬ 
ing effects, produced by a different inclination between the inner 
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velocity (LSR) [km/s] 



2.160 2.162 2.164 2.166 2.168 2.170 2.172 

X [urn] 

Fig. 9. Comparison of our interferometric observations with the inter¬ 
ferometric observables derived from our best disc-wind model P5 plus 
asymmetric disc model. From top to bottom, observed Bry line profile 
(grey) and model line profile (pink), observed visibilities (grey dots with 
error bars) and model visibilities (coloured lines), observed and model 
wavelength-differential phases, and observed and model closured phase. 


discs produced by tidal in t eractions (see e.g.lTuthill et al.ll200ll : 

iDemidova & Grininll2014l: iMarino et al.ll2015lb 

To investigate the origin of the asymmetries in the HD 98922 
disc intensity distribution, we should analyse large- and small- 
scale asymmetries. 

Asymmetries on small scales include i) asymmetric inner 
brightness distribution, due to shadowing caused by the verti¬ 
cally extended structure of the inner rims (which can appear un¬ 
usually asymmetric if the dis c geometry is not exactly face-on; 
see e.g. iMonnier et al.ll2006l) . or ii) brightness peaks and gaps 
along the disc azimuthal angle. 

The large scale asymmetry is represented by the different 
brightness distribution between the NNE and SSW side of the 
disc (whose ratio is about 2.5). This feature is detected at both 
VLTI (AMBER and PIONIER) and SINEONI spatial scales. 

The different emission between forward and backward scat¬ 
tering could be the simplest explanation. The fact that the bright¬ 
ness asymmetry lies approximately in the direction of the system 
axis also supports this hypothesis. On the other hand, this asym¬ 
metry is observed both at large (SINEONI) and small (VLTI) 
scales, and the SINEONI data are likely dominated by scattered 
light, whereas the VLTI data are possibly dominated by thermal 
emission. These considerations are in favour of an intrinsic per¬ 
turbation of the disc. 

In pr inciple, disc warping N the most plausible explanation. 
Recentlv IPemidova & GrininI (l2014l) modell ed the disc asymme - 
try of LkHa 101 (a Herbig Ae star; see e.g. iTuthill et alj|2001h . 
which shows a warped disc. Such asymmetry is very similar 
to tha t observed in HD98922 (see e.g. Eigure 1 in iTuthill et al.l 
1200Ih . The warp is generated by the tidal effects of a low-mass 
companion (with mass ratio q ranging from 0.1 to 0.01) with 
its orbit slightly inclined with respect to the disc plane. The 
warp height varies as the companion extracts matter out of the 
mid-plane of the disc. These perturbations generate an azimuthal 
asymmetry of the extinction. As a result, the illumination of the 
disc by the star becomes asymmetric. 

In the case of HD 98922, the SINEONI data exclude the 
presence of a companion as close as ~20AU from the cen¬ 
tral source down to ~0.5Mo (see Sect. [H i.e. q >0.1). More¬ 
over, the VLTI/PIONIER reconstructed image (see lKluska et al.l 
l2014l) does not show any evidence of a close binary (up to 
a few AUs from the star), but for these data no sensitivity 
limit is prov ided. Indeed, give n the VLTI/PIONIER sensitivi¬ 
ties (see e.g. lAbsil et al.ll201 ih . the dynamical range of the re¬ 
constructed image should be lower than or equal to 100 (Kluska 
priv. comm). This roughly gives a lower limit of H ~ 10.2 mag 
or Mh ~ 2 mag at th e adopted distance of 440 pc. By adopting 
the ISiess et ^ (l2000l) evolutionary tracks (see Sect. HI, we ob¬ 
tain an upper limit for the mass of any undetected companion 
of ~0.4-0.5Mo. Our reasoning implies that the presence of a 
low-mass object (M* < 0.5 Mq) or a massive protoplanet can¬ 
not be ruled out by both SINEONI or VLTI observations. As a 
consequence the hypothesis of a warped disc, originating from 
a very low-mass companion (i.e. with mass ratio in the range of 
0.1-0.01), still holds. 

7.2. Inner circumstellar structure of HD 98922 


Our interferometric observations show that the Bry line visibil¬ 
ities are higher than those of the continuum, indicating that the 

__ Bry emitting region is less extended (~0.3AU in radius) than 

and the outer disc (see e.g. iMarino et al.lF2015h ; d) arcs can be the continuum emitting region (~0.7AU in radius). The pure 
also produced by dust density variations and stellocentric grain- Bry line visibility of HD 98922 is 0.8 at the longest baseline 
size segregation; e) arcs and rings can also originate from warped (89 m), indicating that the Bry originates from a spatially ex- 
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tended region. Our interferometric observables (Bry line profile, 
visibilities, differential phases, closure phases) are sufficiently 
well reproduced by a disc-wind model with a half opening angle 
of 30°, wind footpoints extending from ~0.1 AU to ~1 AU, and 
M„ ~ 2 x 10 -^ Mo yr-^ 

The wind mainly originates from the inner disc region, 
beyond the corotation radius (Rear), which can be placed at 
~0.06AU from the star (assuming v^f i/n/=50+3kms“* and 
i-2Q°, see Table|3. Notably, the wind model extends slightly 
beyond the dust sublimation radius, located at ~0.7 AU. This is 
required to fit both the low visibilities and the large differential 
phases. Indeed, a more compact disc-wind model would produce 
small or no differential phases and it would not re produce our in¬ 
terferometric data (see lGarcia Lonez et alJl20f^ . 

It is worth mentioning that some of the previous HR 
AMBERA'LTI studies of Herbig AeBe s ti ^s (e.g. MWC 297 
and H D 163296; see IWeigelt et al.l 1201 it iGarcia Lonez et ahl 
I2015L respectively) show similar results to those described 
here, namely that the Bry line is emitted from an extended 
region that can be modelled with a disc wind. However, 
MWC 297 (SpT^BlV) and HD 163296 (SpT^AlV) show disc- 
wind geometries that differ from the one modelled in HD 
98922 (SpT=B9V). In the first case, the disc wind was en¬ 
tirely launched from a region located in the continuum disc, 
i.e. well beyond the inn er radius of the continuum-emitting 
disc dWeigelt et al.ll20ll]l . The disc-wind model of MWC 297 
extends from ~17 to 35 R* with a large half opening angle of 
80°. Moreover, the values of the fi and y parameters, which regu¬ 
late the mass load along the streamlines, are quite small (1 and 2, 
respectively). This indicates that th e mass load is shifted towards 
the outer streamlines of the wind. IWeigelt et alJ (1201 Ih suggest 
that such features might be caused by the intense radiation pres¬ 
sure of the central source. On the other hand, HD 163296 has a 
very compact disc wind (extending ~0.16AU ), confined to the 
inner gaseous disc (iGarcia Lonez et al]l2015h . In this case, u)i 
and u)N radii are 2 and 4R, with a half opening angle of 45°. 
In contrast to MWC 297 a nd HD 98922 HD 1 63296 drives a 
well-collimated jet (see e.g. IWassell et alJl200^ : therefore, the 
smaller u values might be related to this. The modelled disc 
wind of HD 98922 has an intermediate geometry, with the Bry 
line-emitting region located mostly in the inner disc, but ex¬ 
tending slightly beyond the dust sublimation radius. Our disc- 
wind model extends from 3 to 30 R, with a half opening angle 
of 30°. Although spectroscopy indicates a s trong outflow origi - 
nates from its disc (e.g. this work; see also lGradv et alj[l996l) . 
no signature of any collimated jet has been detected so far. On 
the other hand, the fi and y parameters are identical to those of 
HD 163296, suggesting that the mass load is similarly distributed 
along the stream lines, but along a larger portion of the disc. 

Lollowing these considerations, our modelling seems to sug¬ 
gest that different types of outflows/jets can be modelled with 
different disc-wind geometries. Moreover, these disc-wind ge¬ 
ometries seem to display an evolutionary trend, following the 
different spectral types and/or evolutionary stages of the central 
sources. It is also tempting to speculate on some additional cor¬ 
relations between, for example, the Bry size, its location, and 
the stellar type and mass of the central source. If our previous 
speculations are correct, we would then expect that disc wind ge¬ 
ometry, kinematics, and Bry size depend on the central source, 
namely its mass and/or evolutionary stage. Nevertheless, it is 
also clear that the sample so far studied is still too limited to 
draw any firm conclusion. 


8. Conclusions 

In this paper, we analyse the main physical parameters and 
the circumstellar environment of the young Herbig Be star HD 
98922. We present AMBLR high spectral resolution (R=12 000) 
interferometric observations across the Bry line, as well as 
UVLS (high-resolution UVB spectroscopy) and SINLONI- 
AO assisted (NIR integral field spectroscopy) ancillary data. 
To interpret our observations, we also developed a magneto- 
centrifugally driven disc-wind model along with a very simple 
asymmetric continuum disc model. The main results of this work 
are the following; 


- UVLS high-resolution spectrum displays Hi photospheric 
absorption lines from the Balmer series (from H 11 to H/1) 
with superimposed emission originating from circumstellar 
activity. Additional weaker photospheric absorption features 
from atomic lines (He i, Ca ii, Le i, Le ii. Mg i. Mg ii. Si i. Si ii) 
and circumstellar emission lines (Le ii and Ti ii) are also de¬ 
tected in the spectrum. 

- Our analysis of the UVLS spectrum indicates that HD 98922 
is a young (5x10^ yr) Herbig Be star (SpT=B9V), located at 
a distance of 440+“ pc, with Macc=(9±3)xl0“^Moyr''. 

- SINLONI fif-band AO-assisted imaging shows a spatially re¬ 
solved circumstellar disc (~140 AU in diameter) with asym¬ 
metric brightness distribution. In particular the NNL side is 
~2.5 times brighter than the SSW side. Other features at 
smaller spatial scales include peaks and dips along the disc. 

- VLTI/AMBLR high spectral resolution observations 
(R=12000) allow us to study the interferometric observ¬ 
ables (i.e. line profile, visibilities, differential phases, and 
closure phases) in many spectral channels across the Bry 
line. Differential phases are V-shaped and line visibilities 
are slightly asymmetric. 

- Our interferometric observations indicate that the Bry line 
visibilities are higher than those of the continuum, i.e. the 
Bry emitting region is less extended than the continuum 
emitting region. Lor the first time, HD 98922 continuum- 
corrected Bry pure line visibilities are 0.8 at the longest 
baseline (89 m), indicating that the Bry originates from 
a spatially extended region. By fitting geometric Gaus¬ 
sian and ring models to the derived pure line visibilities, 
we infer a HWHM Gaussian radius of 0.65+0.09 mas (or 
~0.29+0.04AU at a distance of 440 pc) and ring-fit radius 
of 0.70+0.09 mas (~0.31+0.04AU) for the Bry emitting re¬ 
gion. The Bry line radius is smaller than the inner rim radius 
of 0.7 AU of our continuum temperature model. 

- To obtain a more physical interpretation of our Bry 
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proximately reproduces all the interferometric observables. 
As a result, our modelling suggests that the observed Bry 
line-emitting region mainly originates from a disc wind with 
a half opening angle of 30°, wind footpoints extending from 
~0.1 AU to ~1 AU, and M„,=2xlO'^ Mq yr"'. 

- The observed V-shaped differential phases and slightly 
asymmetric visibilities can be reproduced by combining a 
simple asymmetric continuum disc model with our Bry disc- 
wind model. 
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Appendix A: AMBER interferometric resuits 26 Appendix B: Modeiiing resuits for disc-wind modei 

December 2012 P5 and symmetric disc modei 
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Fig. A.l. AMBER observation of HD 98922 at spectral resolution of 
R=12 000 (26 Dec 2012). From top to bottom: wavelength dependent 
flux, visibilities, wavelength-differential phases, and closure phases ob¬ 
served at three different projected baselines (see labels in figure). For 
clarity, the differential phases of the first and last baselines are shifted 
by -1-30° and -30°, respectively. 


Fig. B.l. Comparison of our interferometric observations with the in¬ 
terferometric observables derived from our best disc-wind model P5 
plus the symmetric disc model described in Sect. 16.21 From top to bot¬ 
tom, observed Bry line profile (grey) and model line profile (pink), 
observed visibilities (grey dots with error bars) and model visibilities 
(coloured lines), observed and model wavelength-differential phases, 
and observed and model closure phase. 


































































































































































